Drude Weight at Finite Temperatures for Some Non-Integrable Quantum Systems in 

One Dimension 
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Using conformal perturbation theory, we show that for some classes of the one-dimensional 
quantum liquids that possess the Luttinger liquid fixed point in the low energy limit, the Drude 
weight at finite temperatures is non- vanishing, even when the system is non-integrable and the total 
current is not conserved. We also obtain the asymptotically exact low-temperature formula of the 
Drude weight for Heisenberg XXZ spin chains, which agrees quite well with recent numerical data. 



Quantum one-dimensional (ID) systems show anoma- 
lous transport properties quite different from higher di- 
mensional systems [1-16]. In particular, Zotos et at 
showed the remarkable fact that in integrable ID sys- 
tems the Drude weight is non-vanishing at finite tem- 
peratures T, and the transport is ballistic, because of 
the presence of non-trivial conservation laws [1,2]. The 
Drude weight is a coefficient of the singular part of the 
conductivity: cr(u>) = ttD(T)S(u>) + (regular part). Their 
proposal has been followed up by extensive studies based 
on various techniques, such as the Bethe ansatz [3-5,9], 
numerical methods [6-8,10], and bosonization [12]. Some 
authors also investigated how this anomalous transport 
is affected by the perturbations that break the integrabil- 
ity [6-8,10]. Recent numerical studies suggest that some 
classes of non-integrable ID systems may also show the 
non-zero Drude weight at finite temperatures, implying 
that the integrability is not a necessary condition for 
D(T) ^ [6-8]. However, the relation between the fi- 
nite Drude weight and non-integrable perturbations has 
not yet been fully elucidated so far. The main purpose of 
this paper is to address this issue. Using conformal per- 
turbation theory, we explore the effects of integrability- 
breaking perturbations on the finite-temperature Drude 
weight in ID Luttinger liquids including spin systems 
and one-component fermion systems. As a by-product, 
we also obtain the asymptotically exact low-temperature 
formula of the Drude weight for Heisenberg XXZ spin 
chains. 

The Hamiltonian consists of the low-energy fixed point 
part Hq, and the irrelevant perturbation H' that may 
render the system non-integrable. Namely, H — Hq+H'. 
Hq is given by the Hamiltonian of the c = 1 Gaussian 
model [17], 



dx 



v[(d x <P(x)) 2 + {d x 0{x)f 



(1) 



charge density (or spin density) operator and the corre- 
sponding current operator, which satisfy the continuity 
equation, are, respectively, given by p(x) = y/Kd x (f>/^/Tr, 
J(x) = — \[Kdt§l \pK. Here K is the Luttinger liquid 
parameter. The irrelevant perturbation H is expressed 
in terms of primary fields of the c — 1 universality class. 
We are concerned with the case that the total current 
J = J dxJ(x) is not conserved by the interaction H': 
[H' , J] ^ 0. This happens when H' contains the cosine 
interaction cos(a<^(x)) which stems from Umklapp pro- 
cesses. In general, H — Hq + H' is not integrable. For 
instance, the multiple-frequency cosine interaction like 
/ dxY^ n [9nCOs(/3 n (j)(x)) + g' n cos(j n 6(x))] breaks the in- 
tegrability of the system [18]. 

To show the presence of non- vanishing Drude weight, 
following Zotos et al. [2], we exploit the Mazur inequality 
which gives the lower bound for the Drude weight. The 
lower bound is expressed in terms of non-trivial conserved 
quantities, which can be found in our case as follows. In- 
troducing a cylindrical geometry with a system size L, 
we write the spatial translation operator as [19], 



I = 



L dx 



T(i)) = ^(L„-I„), (2) 



for the system of the linear size L with the velocity 
v. Here the boson fields <f> and 9 satisfy the canonical 
commutation relation, [4>(x), d x 6{x')] = iS(x — x'). The 



where T(x) (T(x)) is the holomorophic (anti- 
holomorophic) part of the stress tensor. L n and L n 
are the Virasoro generators. For any local field O(x), 
the commutation relation [J, 0(x)] = —id x O(x) holds. 
Note that 0(x) should not be multiplied by a c-number 
function f(x), because I does not operate on f(x) as 
translation [20]. In the case that the Hamiltonian is 
written in terms of local operators, H = J dxO(x), 
we have [I, H] = under the periodic boundary con- 
dition. Thus / is a non-trivial conserved quantity in 
the perturbed system. An important role played by 
the conservation law [I, H] = for transport properties 
was also noticed by Rosch and Andrei before [12,21]. 
For the c = 1 Gaussian model (1), / is expressed as 
I = v J dxd x 9(x)d x 4>(x). This is nothing but the free 
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field representation of the Virasoro generator (2). Us- 
ing the equation of motion of 4>(x) in the presence of 
H', we rewrite I as, I = J dxd t <j>{x)d x <j>{x) + I' with 
I' = —i j dx[H', <j)(x)]d x (j)(x). The total momentum is 
given by P = I + ppJ/K, where pf is the Fermi mo- 
mentum. For the Hamiltonian H = Hq + H', although 
the system is translationally invariant [I, H] = as ex- 
pected for continuum field theory, the charge current 
(or spin current) and the total momentum are not con- 
served, i.e. [P, H] ^ 0, and hence its transport property 
is non-trivial. With the conserved quantity /, the Mazur 
inequality reads [2], 



D(T) > 



1 (JI) 2 
LT (P) 



(3) 



In general, the system has the reflection symmetry 
(particle- hole symmetry in terms of fcrmion fields), (f> — > 
—(j>, which leads to (J I) = 0. However, in the presence 
of a symmetry-breaking field, —hd x (f> (magnetic field in 
spin systems, or chemical potential in fcrmion systems), 
(J I) may be non-zero. We can easily see that this ex- 
ternal field still conserves I. The symmetry-breaking 
field is incorporated into the shift of the boson field 
4>(x) = 4>(x) — hx. Then, we have 



LT LT \K 



with 



-Jr\H; f dujlmIl R (uj)coth(^-) + 
T V K J K J \2TJ 



oj \ {J I') 



LT 



(4) 



n«(o;) = Y J dt([J(q,t),J(-q,0)}) 



e^Vo- (5) 



Here we have used the fluctuation-dissipation theorem in 
the last line of (4). In the case of ImII fl (cij) = 0, the con- 
ductivity is given by Re[er(ii>)] = KvS(lu), showing the 
presence of the finite-temperature Drude weight(see (6) 
below). In the other cases, the first term of the right- 
hand side of (4) is non- vanishing, because Imll^u;) is 
non-negative. Note that the first term of (4) can not be 
canceled with the second term for the following reason. 
Let a the lattice constant, and s — 2 the dimension of 
the leading term of H' . Then, the second term of (4) is 
proportional to a s ~ 2 , of which the order in terms of a 
is different from that of the first term. Hence these two 
terms can not be canceled with each other, and give a 
non-zero lower bound for the Drude weight (3) provided 
that (I 2 )/LT is finite. It is easily seen that (I 2 ) / LT is 
never divergent, since it is given by the derivative of the 
free energy, and the free energy has no thermodynamic 
singularities in ID systems. As a result, we have shown 
that as long as the translational invariance is recovered in 
the scaling limit (i.e. I is conserved), the Drude weight at 
finite temperatures is non-vanishing for the ID systems 



that have the Luttinger liquid fixed point in the presence 
of particle-hole symmetry-breaking fields (magnetic fields 
in spin systems), even if the system is non-integrable, 
and the charge current (or spin current) is not conserved. 
One may question to what extent this argument is rele- 
vant to lattice systems, which are, at short-length scale, 
not translationally invariant. As the low-energy fixed 
point of the systems is the Luttinger liquid, the correla- 
tion length at T = is infinite, characterizing the critical 
state described by the c = 1 universality class. Thus, at 
sufficiently low but finite temperatures, the correlation 
length is much longer than the lattice constant, render- 
ing the systems approximately translationally-invariant, 
provided that there exists no other length scale in the 
Hamiltonian. It should be stressed that, as mentioned 
before, to preserve the conservation law of /, the inter- 
action H' must not contain the c-number function f(x) 
that introduces another length scale [20]. 

Hitherto, our argument is based on the non- 
perturbative analysis. However, it is restricted to the 
case with symmetry-breaking fields. We investigate the 
case without the external fields, by computing the con- 
ductivity pcrturbatively from the Kubo formula, 



a(uj) = 



Kv 



7T 



(6) 



For concreteness, we introduce a specific model given by, 



H'= [dxJ2 



g m cos([3 m 4>(x)). 



(7) 



In general, the multiple-frequency cosine terms (7) break 
the integrabililty of the system [18]. It is noted that the 
ratio of any two of /3 m 's must be rational to ensure that 
the perturbative expansion in term of H' does not gener- 
ate relevant interactions which invalidate the Luttinger 
liquid fixed point. For simplicity, we assume [3 m = m[3\. 
Let us, first, consider the case that g n = for n > 2 and 
only a single cosine interaction exists. In this case, the 
system is the standard Sine-Gordon model, which is in- 
tegrable. We expand the imaginary time current-current 
correlation in term of g\, 



• l/T 



U(iu>) = J dx 

h i dx r dr i d{i) -- An) 
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dT(Td T <j)(x,T)d T (i)(0,0)) e- 

l/T 



(Td T (t>(x,T)d T (t)(0,0)e lai ^e la2 ' p( > 2) ...e la ^ n) ) . (8) 



Here / d(n) = dx n J^ /T dr n , cf>(n) = (j>(x n ,T n ), (■ ■ -) 
is the average with respect to Hq. Using the conformal 
Ward identity, we rewrite the n-th order term of (8) as, 



9i 

2 n nl 



Mi -««)). (9) 
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where 



F aia2 ... an (iuj(l-8 ij )) = J d(l)...d(j-l)d(j + l)...d(n) 

(10) 



G{iuj) 



l dx L' ' 



dT(Td T <j>(x,T)<j>(0,0)) t 



(11) 



with (Td T <p(x,T)<p(0,0)) = i(T/4)coth[(irT/v)(x + 
ivr)] + c.c. For lo 0, G{ito) = v/(iLo), and for lo = 0, 
G(0) = v/T. Correlators of vertex operators exp(mi^) 
in (10) arc non-vanishing only for YJ™ on = because of 
U(l) symmetry. Then, without loss of generality, we can 
put tt2m-i = Pi, CK2m = —01, and n = 2m. It is difficult 
to obtain the explicit expression of (10) as a function of 
lo. However, to show the presence of the Drude weight, it 
is sufficient to investigate the low-frequency behavior of 
F aia2 ... an — Sij)). For this purpose, using the con- 
formal transformation z = cxp(2irTw/v) (w = x + ivr), 
we change the geometry of the system from strip with a 
width 1/T to two-dimensional plane. The correlator of 
vertex operators is transformed as, 

J d{l)...d(n- i)( e ^i^(«'i.»i)... e -i/3i^(«'n,«)„)^ oe -iw(T i -r n ) 
= J d 2 z 1 \z 1 \£~ 2 ...d 2 z n _ 1 \z n _ 1 fi- 2 \z/^e-^-^ 



2irT 



J-n-4n+2 



^ e »/3l^(zi,Zl)___ e -*/3l4>(z n ,«n)^ _ ^2) 



Here we have relabeled the index as n — > j. Following 
Konik and Leclair [22] , we map the IR region to the ultra- 
violet (UV) region by the conformal transformation z' = 
l/z. Note that under this transformation the expression 
of the right-hand side of (12) is unchanged. Thus, al- 
though the IR region of the w-coordinate is mapped to 
both the IR and UV region of the z-coordinate, we need 
only to consider the UV behavior in the z-coordinatc. 
Exploiting the short distance expansion of operator prod- 
uct for the U(l) Gaussian theory [19], we find that the 
dimension of (exp(if3i<p(zi, z\))... exp(—if3i4>(z n , z n )))o m 
the UV region is f3 2 n/(2ir), and hence the dimension of 
the right-hand side of (12) is zero. Therefore, the IR 
singularity of (10) is at most logarithmic, if it exists. 
We can apply the same argument to d k F aia2 ,,, an (iu(l — 
Sij)) I 'd(iLo) k , and find that its dimension in the UV re- 
gion of the z'-coordinate is also zero, because the deriva- 
tive with respect to u does not give rise to an ex- 
tra dimension in the right-hand side of (12). Hence, 



that the expansion of H R (ui) in terms of gi does not give 
rise to singularity stronger than 1/uj. Thus, the Drude 
weight term of a(w) is not eliminated by resummation 
of higher order singularities. A careful treatment of the 
limit lo —* for the second order term of (9) verifies that 
the strongest singularity of <j(uj) is not of order <~ 1/lo 2 , 
but <~ 1/u). Actually, the presence of the double pole 
1/lo 2 is forbidden, because it breaks the non-negativity 
of Recr(oj). Consequently, we have the Drude weight at 
finite temperatures as expected from the integrability of 
the model. 

It is straightforward to generalize the above argument 
to the non-integrable case with multiple cosine interac- 
tions (7). The perturbative expansion gives the same 
expressions as (8) and (9), except that on takes the val- 
ues of /3i, /?2 = 2/3i, /?3 = 3/3i...and so forth. Applying 
the above argument to this case, we also find that the di- 
mension of F aia2 ,,, an (iLo(\ — S^)) and its derivative with 
respect to to is zero, and that F aiCe2 „. an (iu>(l — Sij)) has 
no IR singularity for lj — > 0. This observation leads to 
the conclusion that the 1/w-singularity of a(co) can not 
be eliminated by the resummation of higher-order singu- 
larities. Therefore, the finite-temperature Drude weight 
exists for the non-integrable model (1) and (7). The above 
analysis generalizes Zotos et aUs original proposal par- 
tially to non-integrable systems. This remarkable result 
is physically understood as follows. The existence of non- 
trivial conservation law permits macroscopic numbers of 
level crossing [23]. It is expected that, even for non- 
integrable systems, level crossing may be possible pro- 
vided that a certain non-trivial conservation law holds. 
Then, the presence of macroscopic numbers of degener- 
ate levels results in the non-vanishing Drude weight at 
finite temperatures [6,11]. 
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FIG. 1. Plot of D(T) vs T for f3 2 = 1, 1/12, 1/8, 5/6, 
3/5 (from bottom to top). For comparison, the QMC 
data for [3 2 — 5/6 is quoted from ref. [10] (black circles). 

On the basis of the above perturbative analysis, we cal- 
culate the Drude weight at low temperatures for Heisen- 
berg XXZ chains, 



QlQ2 



.. Qn (w(l - S^)) has no singularity in the limit of H X xz = jy^[S x S x + S v S v + AS Z S Z ] (13) 



lo — > 0. After the analytic continuation in the upper 
half plane iio — > lo + iO, we can expand (10) for small lo: 
F aia2 ... an {(o(l-Sij)) = a +a 1 to(l-5 i j) + a 2 uj 2 {l-5 i: j)--. 
Substituting this into (9), and from YJ i on = 0, we find 



where A = — cos(ir f3 2 ). We consider only the massless 
regime — 1 < A < 1. The low energy effective Hamilto- 
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nian for this regime is completely obtained by Lukyanov 
[24]. Using his result, we obtain the leading temperature 
dependence of the Drude weight at low temperatures. 
For 2/3 < f3 2 < 1, up to logarithmic corrections, it is 
given by, 



D(T) = D(0) 



A = B 



2 



sin(27r//3 2 ) 

87T/3 2 



2tt 2 cot 2 



A 2 (2a7rT)' 



4r(r 2 ) 
r(i-/3- 2 ) 



r(l + /3 2 /(2-2/3 2 )) 



20rF(l + l/(2-2/3 2 )) 



(14) 

(15) 
(16) 



where a = 2(1 - /? 2 )/( Jsin(7r/3 2 )), £>(0) = l/(27ra/? 2 ), 
and B(x,y) = r(x)r(y)/r(a; + y), ^(x) = T'{x)/Y(x). 
For < (3 2 < 2/3, we have, 



D(T) = D(0) 1 



6A_ + A+ 
12 



(a7rT) : 



A + = — tan 

27T 



A_ = 



/? 2 



2 - 2/3 2 / 
r(3/(2-2/3 2 ))r 3 (/3 2 /(2-2/3 2 )) 



12tt T(3/3 2 /(2 - 2/3 2 ))r3(l/(2 - 2/3 2 )) ' 



(17) 
(18) 
(19) 



At the isotropic point /3 2 = 1, logarithmic corrections 
caused by marginal interaction are incorporated into the 
renormalization of the running coupling constant, and we 
end up with, 



D(T) = J -{1 



9 



log(ff) 



+ 0( 5 2 )), 

log(V27 e T +1 / 4 JIT). 



(20) 
(21) 



Here 7 is the Euler's constant. D(T) is non-vanishing 
at the isotropic point in accordance with ref. [10]. The 
above formulae are applicable only in the low temper- 
ature region T <C a -1 . We show the plot of D(T) as a 
function of T for several values of 1 in FIG.l. The result 
for (3 2 =5/6 agrees well with the recent quantum Monte 
Carlo (QMC) data obtained by Alvarez and Gros [10]. 
As 1 decreases, the temperature dependence becomes 
stronger. However, it should be cautioned that, for small 
[3 2 , the applicable temperature range becomes narrower; 
e.g. for [3 2 = 1/12, a -1 = 0.1411. In principle, we can 
improve the formulae by taking into account higher-order 
corrections. 

In summary, we have shown that in a wide class of 
ID integrable and non-integrable systems, as long as the 
low-energy fixed point is the Luttinger liquid, and the 
systems recover translational invariance in the scaling 
limit, the Drude weight is non- vanishing at finite tem- 
peratures, leading to the ballistic transport, even when 
the charge current (or spin current) is not conserved. We 



have also obtained the low-temperature formula of the 
Drude weight for the Heisenberg XXZ chains in the mass- 
less regime. 
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